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<N ! Abstract 

CN ■ We calculate the contributions of the Kaluza-Klein (KK) modes to the Kl — K$ mass 

difference AMk, the parameter ek, the B® s — B® s mixing mass differences AMd lS and rare 
JZL|| decays K + — > tt + vv, Kl — ► n°vv, Kl —>■ /u + /x~, B — > X s ^vv and B s ^ —>■ fifi in the Appelquist, 

On! Cheng and Dobrescu (ACD) model with one universal extra dimension. For the compactification 

scale 1/R = 200 GeV the KK effects in these processes are governed by a 17% enhancement 
of the AF = 2 box diagram function S(xt, 1/R) and by a 37% enhancement of the Z° penguin 
diagram function C(xt/1/R) relative to their Standard Model (SM) values. This implies the 
suppressions of \Vtd\ by 8%, of fj by 11% and of the angle 7 in the unitarity triangle by 10°. 
AM S is increased by 17%. AMk is essentially uneffected. All branching ratios considered in 
this paper are increased with a hierarchical structure of enhancements: K + — » k^vv (16%), 
K L -> ir°vv (17%), B -► A>P (22%), {K L -► /^) SD (38%), £ -> A>z^ (44%), 5 d -> M/2 (46%) 
and S s — ► fj,p, (72%). For 1/R = 250 (300) GeV all these effects are decreased roughly by a factor 
of 1.5 (2.0). We emphasize that the GIM mechanism assures the convergence of the sum over 
the KK modes in the case of Z° penguin diagrams and we give the relevant Feynman rules for 
the five dimensional ACD model. We also emphasize that a consistent calculation of branching 
ratios has to take into account the modifications in the values of the CKM parameters. As a 
byproduct we confirm the dominant O^GpmfR 2 ) correction from the KK modes to the Z°bb 
vertex calculated recently in the large mj limit. 



1 Introduction 



During the last years there has been an increased interest in models with extra dimensions. 
Among them a special role play the ones with universal extra dimensions (UED). In these 
models all the Standard Model (SM) fields are allowed to propagate in all available dimensions. 
Above the compactification scale 1/R a given UED model becomes a higher dimensional field 
theory whose equivalent description in four dimensions consists of the SM fields, the towers of 
their Kaluza-Klein (KK) partners and additional towers of KK modes that do not correspond 
to any field in the SM. The simplest model of this type is the Appelquist, Cheng and Dobrescu 
(ACD) model [1] with one extra universal dimension. In this model the only additional free 
parameter relative to the SM is the compactification scale 1/R. Thus all the masses of the 
KK particles and their interactions among themselves and with the SM particles are described 
in terms of 1/R and the parameters of the SM. This economy in new parameters should be 
contrasted with supersymmetric theories and models with an extended Higgs sector. 

A very important property of the ACD model is the conservation of KK parity that implies 
the absence of tree level KK contributions to low energy processes taking place at scales fi <C 1/R. 
In this context the flavour changing neutral current (FCNC) processes like particle-antiparticle 
mixing and rare K and B decays are of particular interest. As these processes appearing in the 
SM first at one-loop are strongly suppressed, the one-loop contributions from the KK modes to 
them could in principle be important. 

The effects of the KK modes on various processes of interest have been investigated in a 
number of papers. In [1] their impact on the precision electroweak observables assuming a light 
Higgs (run < 250 GeV) led to the lower bound 1/R > 300 GeV. Subsequent analyses of the 
decay B — > X s ^j [2] and of the anomalous magnetic moment [3] have shown the consistency of 
the ACD model with the data for 1/R> 300 GeV. The scale of 1/R as low as 300 GeV would 
lead to an exciting phenomenology in the next generation of colliders [4, 5, 6, 7]. Moreover the 
cosmic relic density of the lightest KK particle as a dark matter candidate turned out to be of 
the right order of magnitude [8]. The related experimental signatures have been investigated in 
ref. [9]. 

Very recently Appelquist and Yee [10] have extended the analysis of [1] by considering a heavy 
Higgs (mu > 250 GeV). It turns out that in this case the lower bound on 1/R can be decreased 
to 250 GeV, implying larger KK contributions to various low energy processes, in particular to 
the FCNC processes. Among the latter only the decay B — > X s ~f has been investigated within 
the ACD model so far [2] and it is desirable to consider other FCNC processes. 

In the present paper we calculate for the first time the B® s — B® s mixing mass differences 
AM,j iS , the Kl — Ks mass difference, the CP violation parameter sk and the branching ratios 
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for the rare decays K + — > 7r + i/p, — ► vr zvi/, jf L — > B — > X Sj rfZ/P and f? S)( 2 — > 

in the ACD model with one universal extra dimension. In the forthcoming paper [11] we will 
analyze the decays B — > X,7, 5 — > X s l + l~ and — > 7r°e + e~. In order to be more general 
we will include the results for 1/R = 200 GeV that is only slightly below the lowest value of 
1/R = 250 GeV allowed by the electroweak precision data. 

As our analysis shows, the ACD model with one extra dimension has a number of interesting 
properties from the point of view of FCNC processes discussed here. These are: 

• GIM mechanism [12] that improves significantly the convergence of the sum over the KK 
modes corresponding to the top quark, removing simultaneously to an excellent accuracy 
the contributions of the KK modes corresponding to lighter quarks and leptons. This 
feature removes the sensitivity of the calculated branching ratios to the scale M s 3> 1 / R 
at which the higher dimensional theory becomes non-perturbative and at which the towers 
of the KK particles must be cut off in an appropriate way. This should be contrasted with 
models with fermions localized on the brane, in which the KK parity is not conserved and 
the sum over the KK modes diverges. In these models the results are sensitive to M s and 
the KK effects in AM S ( j are significantly larger [13] than found here. 

• The low energy effective Hamiltonians are governed by local operators already present 
in the SM. As flavour violation and CP violation in this model is entirely governed by 
the CKM matrix, the ACD model belongs to the class of the so-called models with min- 
imal flavour violation (MFV) as defined in [14]. This has automatically two important 
consequences. 

• The impact of the KK modes on the processes discussed here amounts to the modification 
of the Inami-Lim one-loop functions [15]. This is the function S [16] in the case of AMd )S 
and of the parameter ek and the functions X and Y [17] in the case of the rare decays 
considered. In the ACD model these three functions depend only on m t and the single 
new parameter, the compactification radius R. 

• The unitarity triangle constructed from |Kt&/^c&|> AM^j AM S and the sin 2(3 extracted 
from the CP asymmetry cl^k s ls common to the SM model and the ACD model. That 
is, the -R-dependence drops out in this construction. Which of these two models, if any, is 
consistent with the data can only be found out by analyzing AM^ and AM S separately, 
£k and in particular the branching ratios for rare K and B decays that depend explicitly 
on R. 

Our paper is organized as follows. In section 2, we summarize those ingredients of the ACD 
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model that are relevant for our analysis. In particular, we give in appendix A the set of the 
relevant Feynman rules in the ACD model that have not been given so far in the literature. In 
section 3, we calculate the KK contributions to the box diagram function S and we discuss the 
implications of these contributions for AMk, AM^, AM s , ek and the unitarity triangle. In 
section 4, we calculate the corresponding corrections to the functions X and Y that receive the 
dominant contribution from Z°-penguins and we analyze the implications of these corrections 
for the rare decays K + — > — ► 7r°^z/, — > n + fi~, B — > X s ^yv and -B S)C ; — ► /x/x. In 
section 5, we summarize our results and give a brief outlook. 

Very recently an analysis of AM^ S in the ACD model has been presented in [18]. In the first 
version of this paper the result for the function S found by these authors differed significantly 
from our result with the effect of the KK modes being by roughly a factor of two larger than 
what we find. After the first appearence of our paper the authors of [18] identified errors in 
their calculation and confirmed our result for S. However, we disagree with their claim that the 
reduction of the error on the parameter yj Bs d FB d by a factor of three will necessarily increase 
the lowest allowed value of the compactification scale 1/R to 740 GeV. We will address this 
point at the end of section 3. 



2 The Five Dimensional ACD Model 

The five dimensional UED model introduced by Appelquist, Cheng and Dobrescu (ACD) in [1] 
uses orbifold compactification to produce chiral fermions in 4 dimensions. This is not the case 
in the models described in [19], where all fermions are localized on the 4 dimensional brane. 
However, there are many similarities between these two classes of models, and some of the issues 
discussed in this section have already been presented in detail in [19]. We assume vanishing 
boundary kinetic terms at the cut off scale. We also rely on [1, 20, 21]. 

2.1 Kaluza-Klein mode expansion 

The topology of the fifth dimension is the orbifold S 1 /Z2, and the coordinate y = x 5 runs from 
to 2ttR, where R is the compactification radius. The orbifold has two fixed points, y = and 
y = ttR. The boundary conditions given at these fixed points determine the Kaluza-Klein (KK) 
mode expansion of the fields. 

A scalar field 4> has to be either even or odd under the transformation P5 : y — > —y, and 
therefore 



d^(f) + = for even fields 
4>~ = for odd fields 



► &ty = 0,irR. (2.1) 
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These are von Neumann and Dirichlet boundary conditions respectively at the fixed points. The 
associated KK expansions are 



(f) + (x,y) 



(f> (x,y) 



1 



> 7n s(x) H — ; 6f s(x) cos — 



n=l 

oo 



1 sr a,- t \ ■ n y 

= > 0/ six sin — — , 

n=l 



(2.2) 



where x = x M , fi = 0,1,2,3, denotes the four non-compact space-time coordinates. The fields 
(f)^(x) are called Kaluza-Klein modes. 

A vector field A M in 5 dimensions has five components, M = 0,1,2,3,5. The orbifold 
compactification forces the first four components to be even under P5, while the fifth component 
is odd: 

d b A^ = 



A 5 = 



at y = 0, ttR. 



(2.3) 



Hence, the KK expansion of a vector field is 
1 



A»{x,y) = 



A 5 (x,y) 



V^R >) 



1/ six I cos -— , 

n=l 



b sin — . 

(n)V ' R 

n=l 



(2.4) 



A Dirac spinor ip in 5 dimensions is a four component object. Using the chirality projectors 
Pr/l = (1± 7s)/2, a spinor ^ = (Pr + Pl)^> = ipR + tpL has to satisfy either 

<9 5 ^ = I _ <9 5 ^ 







V>+ = o 



at y = 0, 7rP 



or 



V^ = o 



> at y = 0, ttR. 



(2.5) 



The respective KK mode expansions are 



, , , . 1 . . 1 v-^ / / / \ n y / / \ • n v\ 

(0) + 7^1 ^ v R(n) ^ cos p" + ^ L(n) ^ sm "r J ' 

* v rt=l 

1 -,00 

= (0) + T^R ^ V L(ri) ^ C ° S ~R + ^ R(n) ^ Sm ~R ) 



(2.6) 



The zero-mode is either right-handed or left-handed. The non-zero-modes come in chiral 
pairs. This chirality structure is a natural consequence of the orbifold boundary conditions. 

We can derive Feynman rules for the KK modes by explicitly integrating over the fifth 
dimension in the action: 

2ttR 

S = J d 4 x J dy£ 5 = J d 4 x£ 4 . (2.7) 
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Using the KK mode expansions in (2.2), (2.4) and (2.6), the five dimensional Lagrangian £5 
reduces to the four dimensional Lagrangian £4 which contains all KK modes. The field content 
is arranged such that the zero-modes are the 4 dimensional SM particles, whereas the higher 
modes constitute their KK excitations. Moreover there are additional KK modes that do not 
correspond to any field in the SM. 

2.2 Universal Extra Dimensions 

In the UED scenarios, all fields present in the Standard Model live in the extra dimensions, i.e. 
they are functions of all space-time coordinates. 

For bosonic fields, one simply replaces all derivatives and fields in the SM Lagrangian by 
their 5 dimensional counterparts. There are the U(l)y gauge field B and the SU(2)l gauge 
fields A a , as well as the SU(3)c QCD gauge fields. The Higgs doublet is 




X 2 + iX 1 
tp - ix 3 




(2.8) 



where 



x± = 7! [x± ^ ix " ] • (2 ' 9) 

It is chosen to be even under P5 so it possesses a zero-mode. By assigning a vacuum 
expectation value to that zero-mode with the substitution tp — ► v + H, we can give masses to 
the fermions. Note that we label all parameters of the 5 dimensional Lagrangian with a caret. 
It is convenient to introduce 4 dimensional parameters that are related to their 5 dimensional 



counterparts by factors of \/2kR, see appendix A. 

Fermions living in five dimensions are more involved. In order to write down a Lorentz 
invariant Lagrangian, we need one T-matrix for each space-time dimension to satisfy the Clifford 
algebra 

{T M ,T N } = 2g MN , M,N = 0,1,2,3,5. (2.10) 

The metric g MN = diag(l, —1,-1,— 1,-1) is the natural extension of the flat Minkowskian 
metric for one extra space dimension. For the T-matrices, we take 

P--V, , = 0,1,2,3 
r 5 = i-rs, 

with 75 = Z7°7 1 7 2 7 3 . The 5 dimensional kinetic Lagrangian for spinor fields can now be written 
as [20] 

£5 = i$r M d M ip = ? (i? - 75^ 5 ) ip. (2.12) 
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Integrating over the fifth dimension yields 

n=0 

where the sign of the KK mass term n/R depends on the choice in (2.5). 

Note that the zero-mode remains massless unless we apply the Higgs mechanism. Note 
also that all fields in the 4 dimensional Lagrangian receive the KK mass n/R on account of 
the derivative operator d§ acting on them. These tree-level masses are shifted by radiative 
corrections due to gauge interactions and boundary terms localized at the fixed points [22]. 
Since these corrections are a two-loop effect on the processes considered, we will use the tree- 
level mass relations in our calculations. 

2.3 Gauge fixing and Goldstone mixing 

In the 5 dimensional ACD model, we can use the same gauge fixing procedure as in models in 
which the fermions are localized on the 4 dimensional subspace. We adopt the various gauge 
fixing functional given in [19] and adapt them to the case of the electroweak gauge group 
U(1) Y x SU(2) L with one Higgs doublet: 



g B [B»,B 5 , X 3 } = ^= 
g%[A a »,Alx a ] = -j= 



<v^-e(-f*x a + W) 



(2.14) 



where g' and <?2 are the respective 5 dimensional U(l)y and SU(2)l gauge coupling constants. 
This is a natural extension of the 4 dimensional it^-gauge fixing. 

With the gauge fixed, we can diagonalize the kinetic terms of the bosons and finally derive 
expressions for the propagators. Compared to the SM, there are the additional KK mass terms. 
As they are common to all fields, their contribution to the gauge boson mass matrix is propor- 
tional to the unity matrix. As a consequence, the electroweak mixing angle 6 W is the same for 
all KK modes, and we have 

/ cob0 w -sm8 w \ ( A 3M \ 
~ \ sm9 w cos9 w J \ B M J ' 

W±M = 1 ^ iA2M ^ ^ (2 lg) 

v2 




where 



= sm8 w = - / and c w = cos6 w = - g2 . (2.17) 

V 92 + 9 12 V 92 + 9 12 
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For the /x-components, (2.15) and (2.16) already give the mass eigenstates 



(n) 

A* 

(n) 



,3^ 



s 



(2.18) 



W V2 



The zero-modes and have the masses 



M z 



92+9 



,12 



and 



M W = -j92- 



(2.19) 



The components Z 5 ^ and W^ n ) m i x with the Higgs modes xf n )i while A 5 ^ as defined in 
(2.15) is a mass eigenstate. 

Because of the KK contribution to the mass matrix, the Higgs components an d X 3 with 
n / no longer play the role of Goldstone bosons. Instead, they mix with W§ and Z^ to 
form, in addition to the Goldstone modes and Gf n ^, three additional physical scalar modes 



a" n ^ and a^ n y The former allow the gauge bosons to acquire masses without breaking gauge 
invariance. The propagator terms are diagonalized by the orthogonal transformations 



(n)' 



^0 



G fn) 



1 



M 



Z{n) 



MzX(n) ~ ft Z 5(n) 



«(n) 



M 



Z(n) 



n 



-XU + M z Z 5(n) 



Mi 



W(n) 



± 



W, 



± 



(n) ^> 5(n) 



?! 



>) " Mi 



W(ra) 



- X f n) + Mw w* n) 



± 



(2.20) 
(2.21) 



with M z < n \ and M W ( n ) given in (A. 2). 

For the zero-modes, we can identify X( ) an d X(o) as * ne Goldstone bosons that give masses to 
Z^ and With increasing n, the contributions of Z 5(n ) and W 7 ^) dominate the Goldstone 

modes, while xf n ) an d xf n ) provide the main fraction of and a^ n y The physical Higgs 
does not mix with A 5 ^ n y The latter constitutes an additional unphysical scalar mode which 

turns out to be the Goldstone mode for Af , for n > 1. 

(n) - 

Note that the fields a , a ± and yls do not have zero modes. Consequently, the photon Af. 



(0) 



remains massless. 



2.4 Fermion-Higgs coupling 



The Yukawa coupling of the Higgs doublet to the quark fields is a pivotal part of the Lagrangian 
concerning chirality. Analogous to the SM, we write 



£ qH (x,y) = Q'\vV'(f> + &\ u U'i<j 2 (t>* + h.c. 



(2.22) 
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with the three generation gauge eigenstates 



Q! 




K 
\ u 't ) 



v' 



v. 



(2.23) 



(2.24) 



The SU(2) doublets Q! are odd under P5, while the singlets U' and V are even. Due to this 
assignment, the zero-modes have the same chirality as the quark fields in the SM. 

The fermions receive masses both through spontaneous symmetry breaking and the KK 
expansion as described in section 2.2. In order to diagonalize the Yukawa couplings between 
fermions of equal charge, we apply the biunitary transformation 



Qu = Su Qu, 
W = T u U", 



Q'v = &v Qvi 
V = T v V", 



(2.25) 
(2.26) 



so that the resulting mass matrices 



M v = -^= S^XvTv 
V2 



(2.27) 



are diagonal and their eigenvalues (m u ,m c ,mt) and (m^, m s , nib) are non-negative. This step is 
analogous to the SM and leads to the CKM matrix Vckm = S^Sv which is unique for all KK 
levels. Next we apply to each flavour f = u,c, t, d, s, b the unitary transformation 



U 'kn) 

O" 



f(n) 



-7 5 cosQ /(n) sina /(n) \ / U 
7 5 sina /(n) cosa /(n) / \ Q /(n) 



(2.28) 



and the same expression with U replaced by V for the down- type quarks. The fields Qf( n )i Mf(r, 
and £>/(„) are the mass eigenstates in 4 dimensions. The mixing angle is 



tan 2a i 



m f 



for n > 1. 



n/R 

In (2.28) we have used 75 to get both eigenvalues positive and equal to 



(2.29) 



m f(n) 
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(2.30) 



In phenomenological applications we have n/R> 200 GeV and can therefore set all mixing 
angles to zero except for the top quark. 
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The Yukawa coupling of the leptons to the Higgs doublet is similar to that of the quarks. 
The lepton doublets take the form 

T 

) , (2-31) 



a 



C' 



and the singlets are 

( ei\ 



£' 



(2.32) 



Due to their smallness, we can ignore neutrino masses in heavy flavour physics. Hence, the 
neutrino singlets do not couple to the other particles and are therefore omitted. Due to their 
small Yukawa masses, we can set the mixing angles for all KK excitations of the leptons to zero. 
The lepton masses are given by Eq. (2.30). 

The Feynman rules can be found in appendix A. 



3 B% a - B° d s Mixing and e K 
3-1 Bl s - B». Mixing 



The effective Hamiltonian for AB = 2 transitions in the SM [23] can be generalized to the ACD 
model as follows 



n_,AB=2 



G 



^M^{V* b V tq ) 2 m S{x u l/R) 

6/23 



ai 5) (/4,) 
1 + ^^ J5 



Q(AB = 2) + h.c. 



Here /z 6 = 0{m b ), J 5 = 1.627, 

Q(AB = 2) = (bq) V - A {bq)v-A, q = d,s 
with (bq) V - A = &7 M (1 - 75)9 and [23, 24] 

7] B =0.55 ±0.01 

describes the short distance QCD corrections to which we will return below. 
The function S(xt,l/R) is given as follows 



S(x t , l/R) = S (x t ) + ^2s n (x t ,x r , 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



n=l 
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where (rrit = mt(mt)) 



x t 



2 

mf 



~ M 2 



m n = 



n 
R 



and 



So(x t ) 



Ax t - llx 



3xf In xt 



(3.5) 



(3.6) 



4(1 -x t f 2(1 -x t ) 3 
results from the usual box diagrams with (W ± ,t) and (G ± ,t) exchanges with the m t independent 
terms removed by the GIM mechanism. 

The KK contributions are represented by the functions S n (xt,x n ) that are obtained by 
calculating the box diagrams in fig. 1 with a^, G^y Qi(n) an d^4( n ) (i = u,c,t) exchanges 

and multiplying the result by i/4, where 1/4 is a combinatorial factor. We neglect momenta and 
masses of external quarks. Denoting the contribution of the sum of the diagrams corresponding 
to a given pair (m i(n ), m j(n) ) to S n (x t ,x n ) by F(x^ n) ,x j ^ n) ) with 



W ± ,G ± ,a ± 




Ui, QiMi Jk 




W ± ,G ± ,a ± 



W ± ,G ± ,a ± 



W ± ,G ± ,a 



Figure 1: Box diagrams contributing to S n (x tl x n ). We suppress the KK mode number. 



rn 



X i(n) M 2 



i(n) 



(3.7) 



W(n) 



and mj( n ) and M w ^ defined in (2.30) and (A. 2), and using the unitarity of the CKM matrix, 
we have 



S n (x t ,x n ) = F(x t{n) ,x t{n) ) +i ? (x u(n ),x u(n) ) - 2F(x t ( n) ,x u{n) ) . 



(3.8) 



As with increasing n the modes Qt( n ), Ut{n)i Qu(n) an d W u ( n ) become increasingly degenerate in 
mass and 



Xt(n) > x u(n) > 1) 



(3.9) 



the functions S n (xt,x n ) decrease with increasing n so that only a few terms in the sum in (3.4) 
are relevant. This is seen in fig. 2 (a) where we show S n (xt,x n ) as a function of n/R. We will 
discuss this in more detail in section 4.9. 
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The contributions from different sets of diagrams to the functions F(x i r n \,Xjr n \) are collected 
in appendix B. It turns out that the contribution from the pair (ai\,a7^\) is by far dominant. 
We illustrate this in fig. 2 (a). In phenomenological applications it is more useful to work with 
the variables xt and x n than with x^ n y We find 



4(x* - lfxt 



6x n xt — 5xf — 12x n x\ + 15xj + 10x n x 3 — llxj — 4:X n xj + xf 



2x n (x t - l) 3 (3x n + 3x n x t - x t ) In ■ 



l + x n 

6x n x 2 - 2x\ + 14x n x t 3 - 2x\x\ + 6x\ - 2x n xf] In X ' + ' " 



+ ( - 6x 2 n + 2x n x t + 12x^x t (3.10) 



1 + Xr 



In fig. 2 (b) we plot S(xf, 1/R) versus 1/R. For 1/R = 200 GeV we observe a 17% enhancement 
of the function S with respect to its SM value given by So(xt). For 1/R = 250 GeV this 
enhancement decreases to 11% and it is only 4% for 1/R = 400 GeV. 



0.4 
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(a) 



(b) 



Figure 2: (a) Contribution S n of the n KK mode to S(xt, 1/R)- The contributions with a dominate, 
those with only G ± and 1T ± are negligible and not shown, (b) The functions S(xt, 1/R) and So(xt). 



Proceeding as in the SM, we can calculate the mass differences AM q by means of 
Gl 



AM, 



</ " -^VBm Bq (B Bq F 2 Bq )M^S(x t ,l/R)\V tq \ 2 , 



(3.11) 



where Fg is the i? g -meson decay constant and B q the renormalization group invariant parameter 
related to the hadronic matrix element of the operator Q(AB = 2), see [25] for details. This 
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implies two useful formulae 



AM d = 0.50/ps 



B B d F B . 



230 MeV 



\v t 



Id] 



7.8 • IO- 3 



VB 



0.55 



S(x t , l/R) 
2.34 



(3.12) 



and 



AM S = 18.4/ps 



B Ba F Ba 



270 MeV 



\v t 



1 2 



Is 



0.040 



VB 



0.55 



S{x t ,l/R) 



2.34 



(3.13) 



The implications of these results for \V t d\, AM S and the unitarity triangle will be discussed 
below. 

Finally, a few comments regarding the QCD factor t]b are in order. This factor as given in 
(3.3) has been calculated within the SM including NLO QCD corrections that are mandatory 
for the proper matching of the Wilson coefficient of the operator Q(AB = 2) with its hadronic 
matrix element represented by the parameter Bb 3 d and calculated by means of non-perturbative 
methods. As the top quark and the KK modes are integrated out at a single scale = 
0(mt, l/R), the contributions to t)b from scales lower than [i t are the same for the SM and the 
ACD model. They simply describe the finite renormalization of Q(AB = 2) from scales O(nt) 
down to the scales 0(rrib). The difference between QCD corrections to the SM contributions 
and to the KK contributions arises only in the full theory at scales [it = 0(m t ,l/R) as the 
unknown QCD corrections to the box diagrams in fig. 1 can in principle differ from the known 
QCD corrections to the SM box diagrams [23, 24] that have been included in t]b- As the QCD 
coupling constant a s (nt) is small and the QCD corrections to the SM box diagrams of order of a 
few percent, we do not expect that the difference between the QCD corrections to the diagrams 
in fig. 1 and to the SM box diagrams is relevant, in particular in view of the fact that the KK 
contributions amount to at most 17% of the full result. 

For fi » fi t , a s (/i) in the ACD model becomes larger [26] and the QCD corrections to KK 
modes with n » 1 could in principle be substantial. However, as seen in fig. 2 (a), these heavy 
modes give only a tiny contribution to S(x t , l/R) and can be safely neglected. 



3.2 e K 

The effective Hamiltonian for AS* = 2 transitions is given in the ACD model as follows 



H 



AS=2 

eS 



G 



p M^r [XlmSoixc) + >nmS{x u l/R) + 2\ c \ m S Q {x c ,x t )} 

-2/9 



16vr 2 



1 H 'a ^3 



47T 



Q{AS = 2) + h.c. 



(3.14) 
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(3) 

where Aj = V* s V id , a: s is the strong coupling constant in an effective three flavour theory and 
J 3 = 1.895 in the NDR scheme [23]. In (3.14), the relevant operator 

Q(AS = 2) = (sd) v ^ A (sd) V -A, (3.15) 

is multiplied by the corresponding Wilson coefficient function. This function is decomposed into 
a charm-, a top- and a mixed charm-top contribution. The SM function So(x c ,xt) is defined by 

5 (x c , x t ) = F(x c , x t ) + F(x u , x u ) - F(x c , x u ) - F(x t , x u ), (3.16) 

where F(xi,xj) is the true function corresponding to the box diagrams with exchanges. 
One has 



S (x c ,x t ) = x c 



xt 3xt 3xf In xt 



In — - 



(3.17) 



4(1 -x t ) 4(1 -x t ) 2 _ 

where we keep only linear terms in x c <C 1, but of course all orders in xt- 

In view of the comments made after (3.8) and the structure of (3.16), the impact of the 
KK modes on the charm- and mixed charm-top contributions is totally negligible and we take 
into account these modes only in the top contribution that is described by the same function 
S(xt, 1/R) as in the case of AM q . This also means that the Kl — Kg mass difference, AMk, 
being dominated by internal charm contributions in (3.14) is practically uneffected by the KK 
modes. 

Short distance QCD effects are described through the correction factors rji, 772, % and the 
explicitly a s -dependent terms in (3.14). The NLO values of r\i are given as follows [27, 23, 28]: 

771 = 1.45 ± 0.38, 772 = 0.57 ±0.01, r/ 3 = 0.47 ± 0.04 . (3.18) 

The standard procedure allows now to calculate the CP-violating parameter ek [25] 

e K = C £ B K lm\ t {ReA c [rnS {x c ) - rj 3 S (x c ,x t )] - Re\ t 7i 2 S(xt, V-R)} exp(i7r/4) , (3.19) 

where C e = 3.837 • 10 4 is a numerical constant. Bk is the renormalization group invariant 
parameter related to the hadronic matrix element of the operator Q(AS = 2), see [25] for 
details. 



3.3 Unitarity Triangle in the ACD Model 

What is the impact of the KK contributions to the function S on the elements of the CKM 
matrix and in particular on the shape of the unitarity triangle? In order to answer this question 
let us recall a few aspects of the unitarity triangle (UT) shown in fig. 3 and of the Wolfenstein 
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parametrization [29] as generalized to higher orders in A in [30]. The apex of the unitarity 
triangle is given by [30] 



A 2 A 2 
£ = p(l-y), fj = ri(l-—). 

Here A, A, g and r] are the Wolfenstein parameters [29]. Moreover, one has 

V us = X + 0(X 7 ), V ub = AX 3 (g-i V ), V cb = AX 2 + 0(A 8 ), 

Vu = -AX 2 + ^AX 4 [l - 2( e + ir})], V td = AX\l -g-ifj). 



A = (p, f)) 




(3.20) 

(3.21) 
(3.22) 



C = (0,0) 



Figure 3: Unitarity Triangle. 
The lengths R b and Rt are given by 



Rh = 



Rt = 



I Vu d Vub\ 

\y td Vtt\ 
K d v; b \ 



\J g 2 + rj 2 



V 2 ; A 



V ub 



V td 



v cb 



(3.23) 



(3.24) 



and the angles (3 and 7 of the UT are related directly to the complex phases of the CKM elements 
Vtd and V ub , respectively, through 



Vtd = Wdle-^, V ub = \V ub \e~ i \ 
The five constraints on the UT that we have at our disposal at present are: 



(3.25) 



• The R b Constraint: As seen in (3.23), the length of the side AC is determined from 
\V ub /V cb \. 

• Hyperbola (Indirect CP Violation in Ki — > tt-k) obtained from (3.19) and the experi- 
mental value for ek (A = 0.221): 



fj [(1 - g)A 2 V2 S(x t , l/R) + P c (s)] A 2 B K = 0.214 , 



(3.26) 
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where P c (s) = 0.28 ± 0.05 [27, 28] represents the charm contribution that is not affected 
by the KK contributions. 



B% - B^-Mixing Constraint (AM d ): 

\Vtd\ 



A|y, 



cb| 



7.8 • 10" 3 



with 



\V U \ = 7.8 • 10- 



230 MeV 



"0.041" 










/0.55 


0.50/ps 


V VB 



2.34 



5(x t , 1/12) 



(3.27) 



(3.28) 



where AM d = (0.503 ± 0.006)/ps [31]. 
5° - £?°-Mixing Constraint (AM d /AM s 



Rt = 0.87 



AM d 18.4/ps 



0.50/ps V AM S 



1.18 



B dF~B d 



where AM S > 14.4/ps at 95% C.L. [31]. 



The direct measurement of sin 2/3 through the CP asymmetry a(ipKs) in B d 
is not affected by the KK contributions. 



(3.29) 



ipKs that 



The main uncertainties in this analysis originate in the theoretical uncertainties in the non- 
perturbative parameters Bx and y B d FB d and to a lesser extent in £ [32] : 



B K = 0.86 ± 0.15, \/B~ d F Bd = (235^) MeV, £ = 1-18±8;J|. 



(3.30) 



Also the uncertainties in |Vu6/V^fe| and y B s Fb s = (276 ± 38) MeV are substantial. The QCD 
sum rules results for the parameters in question are similar and can be found in [33]. 

With these formulae at hand let us enumerate a few general properties of the values of the 
CKM elements within the ACD model. These are: 

• \Vus\, \Vcb\ an d iKtfcl are usually determined from tree level decays. As in the ACD model 
there are no KK contributions at the tree level, the absolute values of these three CKM 
elements are to an excellent approximation the same as in the SM model. From the point 
of view of the unitarity triangle (UT), this means that the lengths of its two sides, AC and 
CB are common to the SM model and the ACD model. In our numerical analysis we will 
use, as in [34], 

\V ch \ = (40.6 ±0.8) • 10" 3 , 



IV I = A = 0.221± 0.002 



\Vub\ 

\v cb \ 



0.089 ± 0.008, 



\V uh \ = (3.63 ±0.32) • 10" 



(3.31) 
(3.32) 
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• The angle [5 of the UT has been determined recently by the BaBar [35] and Belle [36] 
collaborations from the CP asymmetry a^K s in -B ^ ipKs with a high accuracy giving 
the world grand average [37] 

(sin 2(3)^ Ks = 0.734 ± 0.054 . (3.33) 

As there are no new complex phases in the ACD model beyond the KM phase, the angle 
as extracted by means of cl^k s ls common to both models in question. 

• A similar comment applies to |Vtd| or equivalently the length R t of the side AB in fig. 3, 
when R t is extracted from the ratio AMd/AM s that is independent of S(xt, as seen 
in (3.29). 

Thus when | V^ft/V^bl , sin 2(3 from a$Ks and AM^/ AM S are used to construct the UT, there is 
no difference between the SM and the ACD model as all explicit dependence on 1/R cancels out. 
This universal UT (UUT) [14] that is valid for all MVF models, as defined in [14], has recently 
been determined [34, 38]. 

Now, even though there exists a UUT common to the SM and the ACD model, in view of 
the fact that S(xt, ^ So(xt), only one of these two models, if any, will have £k, AMd and 
AM S that agree with the experimental data. Let us consider AM S first. As seen in (3.22) |Vt s | is 
very close to |V^| because of CKM unitarity. Therefore it is common with an excellent accuracy 
to both models and consequently 

(AM S ) ACD = S(x t ,l/R) 

(AM S ) SM S (x t ) • { ' ' 

However, this ratio is at most 1.17 and the distinction between these two models will only be 
possible provided Bb s F Bs can ' 3e calculated to better than 10% accuracy. A very difficult task. 

The fact that S(xt, 1/R) > So(xt) implies also that with the experimentally known values of 
£k and AMd, one has 

|U td |ACD < \V td \sM, 7ACD < 7SM (3.35) 
as can easily be inferred from (3.26) and (3.28). In particular (3.28) implies 



\Vtd\ ACD I S (x t ) . . 

I^Ism \js(x t MRy 

Thus |Vtrf|ACD can be smaller than |Vtd|sM by at most 8%. In order to determine such a differ- 
ence, more accurate information on the unitarity triangle and the non-perturbative parameters 
entering ek and AM^^ is necessary. Similarly 7acd can be smaller than 7sm by at most 10 
degrees. 
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We illustrate these properties in fig. 4, where we show the 1/R dependence of \Vtd\, 7, Q 
and 77. In obtaining these results we used the following procedure. First using the central 
value |Vt(i|sM = 0.00815 from the SM fit of [34] and rat = 167 GeV we determined (V^Iacd by 
means of (3.36). The result is shown in fig. 4 (a). Using next the central values in (3.31) and 
(3.32) and the 1/R dependence of {Rt) acd = |^td|ACD/(A|^&|), we determined 7, g and fj as 
functions of 1/R. We show in fig. 5 the unitarity triangles corresponding to the ACD model 
with 1/R = 200 GeV and the SM model. 



Figure 5: Unitarity Triangle in the ACD model for 1/R= 200 GeV and in the SM. 

In summary, the CKM elements in the ACD model extracted from AF = 2 transitions and 
the CP asymmetry a^K s are not expected to differ substantially from the corresponding values 
found within the SM. This is very fortunate as the most recent fits of the UT [34, 37, 39, 40, 41, 42] 
based on the SM expressions for AF = 2 transitions agree very well with the direct measurement 
of the angle f5 by means of a^x s - With improved data on a^K s an d in particular AM S and 
improved values for Bk and Bb q F^, a constraint on the compactification radius R from AF = 2 
processes is in principle possible. However, for the time being even the lowest value 1/R = 
200 GeV considered by us is consistent with the present fits of the UT. Setting 1/R = 200 GeV 
and repeating the analysis of [34], that uses the bayesian method [39], one finds [43] the values 
for (g,fj), sin2a, 7, AM S and \Vtd\ m the second column of table 1. For a comparison we give 
the corresponding ranges in the SM. To this end all input parameters of [34] have been used. 

Comparing the two columns in table 1, we observe all the patterns shown in fig. 4. However, 
due to substantial uncertainties in the input parameters, the effects of the KK modes are partly 
washed out. In particular, the suppressions of \Vtd\ and fj amount to 4 — 5% rather than O(10%) 



that we found in fig. 4. This is easy to understand. With no uncertainty in y Bs d FB d , the value 
of \Vtd\ is strongly correlated with S(xt,l/R) by means of (3.12) as AMd is known very well 
experimentally. Once the uncertainties in y BB d FB d are taken into account, the enhancement 



A = {p,fj) 




c=(o,o) 



B = (1,0) 
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of the function S in the ACD model can be compensated by the decrease of both \Vtd\ and 
\J Bb a Fb a , implying a smaller suppression of \Vtd\ than found in fig. 4. Similar comments can 
be made in connection with other quantities in table 1. 

Clearly, low energy non-perturbative parameters as \J Bb a Fb a , Br and £ do not depend on 
the compactification scale. On the other hand they are subject to uncertainties and it is not 
surprising that the fit of the unitarity triangle within the ACD model prefers lower values of 
\J~B~B d FB d and Bk than in the SM. The true effects of the KK modes can then only be clearly 
seen as in fig. 4 when the input parameters have no uncertainties. This exercise shows very 
clearly the importance of the reduction of the theoretical uncertainties in connection with the 
search for new physics. 



Strategy 


ACD (1/R = 200 GeV) 


SM 


fj 


0.342 ± 0.027 


0.357 ± 0.027 




(0.288 - 0.398) 


(0.305 - 0.411) 


P 


0.197 ± 0.047 


0.173 ± 0.046 




(0.102 - 0.296) 


(0.076 - 0.260) 


sin 2a 


-0.23 ±0.25 


-0.09 ±0.25 




(-0.70 - 0.27) 


(-0.54 - 0.40) 


7 


59.5 ± 7.0 


63.5 ± 7.0 


(degrees) 


(45.3 - 74.8) 


(51.0 - 79.0) 


AM S 


is.6±i:§ 


18.0^ 


(ps- 1 ) 


(15.7-26.2) 


(15.4-21.7) 


m (io- 3 ) 


7.80 ± 0.42 


8.15 ± 0.41 




(6.96 - 8.69) 


(7.34 - 8.97) 



Table 1: Values and errors for different quantities in the ACD model with 1/R= 200 GeV [43] and in 
the SM from [34]. The 95% probability regions are given in brackets. 



This discussion makes it clear that it is impossible to claim, as done in [18], that the reduction 
of the error on the parameter \J Bb a Fb a by a factor of three will necessarily increase the lowest 
allowed value of the compactification scale 1/R. The lower bound on 1/R from AM^ depends 
necessarily on the actual value of \J BB d Fs d - With decreasing \J Bb a Fb a the lower bound on 

1/R becomes weaker. One can easily check that decreasing the central value for \J Bb a Fb a 
to 200 MeV, still within the present uncertainties, no significant lower bound on 1/R can be 
obtained even if the error on \J Bb a Fb a is decreased by a factor of three. 

In order to illustrate this point we show in fig. 6 the correlation between 1/R and J Bb a Fb a 
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Figure 6: Correlation between 1/R and \J B Bd FB d for different values of \V t 
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'td\ 



for different values of \Vtd\ and mt = 167 GeV. To this end we have used (3.12) with AM^ = 
0.503/ps. Clearly on the basis of AMd alone it is impossible to say anything about 1/R. However, 



even if \ V t d\ is determined through AM^/AM S , the values of 1/R depend sensitively on y BB d FB d 



with the maximal and minimal values of y Bs d FB d corresponding to maximal and minimal 



values of 1/R, respectively. Thus a significantly improved lower bound on y B-B d FB d could in 
principle provide an improved lower bound on 1/R. However, if future lattice calculations will 



find values for y BB d Fs d that are smaller than the present central value, 1/R could be forced 
to be low in order for the ACD model to fit the value of AMd- 



4 Rare K and B Decays 
4.1 Preliminaries 

We will now move to discuss the semileptonic rare FCNC transitions K + — > -k + vv, K\, — > tt°i/9, 
Kl — ► /U/2, B — > X s ^dyv and B St d — > fJ-fi- Within the SM and the ACD model these decays 
are loop-induced semileptonic FCNC processes governed by Z°-penguin and box diagrams and 
described by two functions X(xt,l/R) and Y(x t ,l/ R) for the decays with vv and up, in the 
final state, respectively. 

A particular and very important virtue of these decays (with the exception of Kl — » /j,p) is 
their clean theoretical character [25] that allows to probe high energy scales of the theory and in 
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particular to measure Vtd and ImAt = lmV t *Vtd from K + — > 7r + i/z^ and i^L — ► 7r°^t> respectively. 
Moreover, the combination of these two decays offers one of the cleanest measurements of sin 2(3 
[44], see [44, 45] for more details. 



4.2 Effective Hamiltonians for K — > nuu and B —*■ X s vu 

The effective Hamiltonian for K + — > 7r + ^P is given in the ACD model as follows 
Gf a 



7~Lefi — 



y/2 2vr sin 2 6 W , ^ 

(=e,/i,T 



^ (K^^l + VT s ^x^(x t , 1/i?)) {sd)v- A {nvi)v-A ■ (4.1) 



The index l=e, fj,, r denotes the lepton flavour. The dependence on the charged lepton mass 
resulting from the box diagram is negligible for the top contribution. In the charm sector this 
is the case only for the electron and the muon but not for the r-lepton. In what follows we will 
set the QCD factor i]x [46, 47, 48] to unity as for m t = mt(m t ) it equals 0.994. 
The function X(x t ,l/R) relevant for the top part is given by 



X(x t , 1/R) = C(x t , l/R) + B up (x t , l/R). 



Here 



C(x t ,l/R) = C (x t ) + J2 C n(xt,X n ) 

n=l 

results from Z°-penguin diagrams with the SM contribution Co(xt) given by 



(4.2) 



(4.3) 



C (x t ) 



x t — 6 3x t + 2 



+ 



In xt 



(4.4) 



x t -i ' (x t -iy 

The sum in (4.3) represents the KK contributions that are calculated from the Feynman diagrams 
in fig. 7 as discussed below. Next 



BT{xt, l/R) = -AB {x t ) + Bf(x t ,x n ) 



(4.5) 



n=l 



results from box diagrams with the SM contribution given by the first term and 



Bo(xt) = ~ A 



x t x t In x t 



(4.6) 



l-x t ' (x t -l) 2 _ 

The sum in (4.5) represents the KK contributions that are calculated from the Feynman diagrams 
in fig. 9 as discussed below. 

The expression corresponding to X(xt,l/R) in the charm sector is the function ^^L- ^ 
results from the NLO calculation [49] and is given explicitly in [48] where further details can 
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be found. As in the case of the charm contributions in the AS* = 2 Hamiltonian, here the KK 
contributions are also negligible. The numerical values for X^ L for \i = m c and several values 
of and m c (m c ) can be found in [48]. For our purposes we need only (A = 0.221) 



Po(X) = 



0.41 ± 0.06 , (4.7) 



where the error results from the variation of A^ti- and m c (m c ). 

MS cv c/ 

In the case of — > t\^vv that is governed by CP-violating contributions only the top 
contribution in (4.1) matters. Similarly the effective Hamiltonians for B — > X s ^vv are obtained 
from (4.1) by neglecting the charm contribution and changing appropriately the CKM factor 
and quark flavours. In all these decays the KK modes contribute universally only through the 
function X(xt,l/ R). 

4.3 Effective Hamiltonians for B s d — > fip, and Kl — ► MM 

The effective Hamiltonian for B s — > l + l~ in the ACD model is given as follows: 

Gt? ol — - 

H eS = --= - . 2fl V t tV ts ri Y Y(x t , l/R)(bs)v- A (U)v-A + h.c. (4.8) 
V 2 2ir sin 

with s replaced by d in the case of Bd — > The charm contributions are fully negligible 

here. In what follows we will set the QCD factor r/y [46, 47, 48] to unity as for m t = rnt(m t ) it 
equals 1.012. 

The function Y(xt,l/R) is given in the ACD model by 
Y(x t , l/R) = C(x t , l/R) + B^(x t , 1/R), (4.9) 
with C(xt,l/R) given in (4.3) and 

oo 

B^(x t , l/R) = -B (x t ) + B^(x t , x n ) (4.10) 

71=1 

with the first term representing the SM contribution. The second term results from the box 
diagrams in fig. 9. 

The effective Hamiltonian for the short distance contribution to Kl — > fifi is given by (4.8) 
with the appropriate change of the CKM factors and quark flavours and a small SM charm 
contribution Ynl [48] in analogy to X^ L in (4.1). The relevant branching ratio will be given 
below. 
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4.4 Z°-Penguin Diagrams 

The function C n (xt,x n ) can be found by calculating the vertex diagrams in fig. 7 and adding 
an electroweak counter term as discussed in detail in [50]. The latter is found by calculating 
the self-energy diagrams of fig. 12 that describe flavour non-diagonal propagation of quark 
fields and subsequently rotating the quark fields appropriately so that this flavour non-diagonal 
propagation does not take place in the new fields. 




Figure 7: Penguin diagrams contributing to C n (xt,x n ). The analytical expressions are listed in appendix 
C. (8) includes the additional diagram with W ± and (a ± ,G ± ) interchanged. 

In the 't Hooft-Feynman gauge for the and propagators, the contribution of the 
diagrams in fig. 7 to the flavour-changing vertex T^ z including the electroweak counterterm is 
given by 

= ' I67/L 9 W V ts V tdCn{x U Xn)s^(l ~ J 5 )d. (4.11) 

We neglect the external momenta in fig. 7 and the masses of external quarks. The function 
C n (xt,x n ) is defined through 

C n (x t , x n ) = F(x t(n) ) - F(x u(n) ) (4.12) 

with the functions F(x t ^) and F(x u ( n )) representing the contributions of the Qt( n ), Ut{ n ) an d 
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Qu(n), U u {n) modes, respectively, 

8 2 

F(x t{n) ) =J2F i (x t{n) ) + (\ - l -s 2 \ J2ASi(x t{n) ). (4.13) 
i=i ^ ' i=i 

Here s w = sin# w , Fj stand for the contributions of diagrams in fig. 7 and ASi denote the 
electroweak counter terms. 

The explicit contributions of various sets of diagrams to these functions are given in appendix 
C. Adding up these contributions we find 



8(x t - 1)2 



x\ - 8x t + 7 + (3 + 3x t + 7x n - x t x n ) In X ' 



(4.14) 



1 + x n 

In fig. 8 (a), we show C n (x t , x n ) as a function of n/R. In this case the convergence of the sum of 
the KK modes is significantly improved by the GIM mechanism so that only a few first terms in 
the sum in (4.3) are relevant. We will return to this at the end of this section. We observe that 
in contrast to the function S of section 3 the diagrams involving only play the dominant 
role among the KK contributions. 

In fig. 8 (b) we plot C{x t , 1/R) versus 1/R. For l/R = 200 GeV we observe a 38% enhance- 
ment of the function C with respect to its SM value given by Co(xt). For 1/R = 250 GeV this 
enhancement decreases to 26% and it is 11% for 1/R = 400 GeV. The significant enhancement 
of C is the origin of the enhancements of the branching ratios discussed below. 

4.5 Z°bb Vertex in the large m t limit 

The result in (4.14) can also be used to find the dominant KK contribution to the Z°bb vertex in 
the large mt limit. As discussed already in [50], in general the calculation of a low energy effective 
flavour violating Z°ds vertex cannot be directly compared with the full calculation of the Z°bb 
vertex. However, as pointed out there in the special limiting case mt S> M\y the two approaches, 
the direct evaluation of on-shell diagrams and the operator product expansion considered by us, 
are equivalent. The reason for this is that in the limit m t S> Mjy with 1/R > mt all the other 
mass scales as m& and external momenta of 6-quarks become negligible compared with mt and 
1/R. Consequently with the definition 

^ = l ^ lr{9LPL + 9RPR)h (415) 

the onedoop contributions to the coupling can be found in the large mt limit by using the 
relation [50] 

fe = ^-fC( I( ,l/JJ). (4.16) 

V2 TT Z 
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In the case of the SM contributions this relation has been already analyzed in detail in [50] 
including a s to the one-loop Z°bb vertex. For the KK contribution we find from (4.14) 

00 / 4 2 2 2\ / 4 4 fi 4 2 2 4\ 

AC= > CJx t , x n ) = — s-g- H — -R + — — n — - — ) R + ... 

^ y ' V96m^ 96 J V4320 2160m^ 864 / 

(4.17) 

with the dominant contribution represented by the first term, see appendix D for the derivation. 
Retaining only this term and using (4.16) we find 

Sg KK = ^ m 4^! (418) 

that agrees with a recent direct calculation of the KK contributiond to the Zbb vertex in [51], see 
their formula (3.6). Interestingly, while the corrections to the asymptotic result (4.18) relevant 
for the Zbb vertex have been found in [51] to be substantial, the result 

AC = ^rR 2 (4.19) 
96M^ v ' 

give a good approximation to the full KK contribution to the function C(xt,l /R) even for 
m t « 167 GeV. 
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4.6 Box Diagram Contributions 



The functions B^ u (xt,x n ) and B n H1 (xt,x n ) can be found by appropriate rescaling of the box 
diagrams contributing to AF = 2 transitions that we considered in section 3. It turns out that 
the box contributions of the KK modes are tiny. For instance for l/R = 200 GeV and m t = 
167 GeV we find Bf(x t ,xi) = 0.0098 and B^(x t ,x{) = 0.0049 with even smaller values for 
n > 1 and larger l/R. Consequently, these contributions can be safely neglected in comparision 
with C n . For completeness we give in appendix D the analytic formulae for B^ u (xt,x n ) and 

B^(x t ,X n ). 



W ± ,G ± ,a ± 



W ± ,G ± ,a ± 





ui, Qi,Ui ^ 



Figure 9: Box diagrams contributing to and to B^f . 



4.7 The Functions X and Y 

Neglecting the box contributions of the KK modes we find 

oo 

X(x t , l/R) = X (x t ) + £ C n (x t , x n ) = X (x t ) + AX 



(4.20) 



n=l 



Y(x t , l/R) = Y (xt) + J2 °n(xu x n ) = Y (x t ) + AF 



(4.21) 



n=l 



with (mt = 167 GeV) 

x t 



X (x t ) 



Yo(x t ) = ? 



x t + 2 3x t - 6 

+ t — hi x t 



x t -i (x t -iy 



x t -4 



+ 



In x t 



1.526 



0.980 



(4.22) 



(4.23) 



[x t -l ' {xt-lf 

summarizing the SM contributions and AX = AY representing the corrections due to KK 
modes. 

In fig. 10 we plot X(xt,l/R) and Y(x t ,l/R) versus l/R. We observe that due to the 
inequality Xq > Yq the relative impact of the KK modes is larger in the function Y. For l/R = 
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(a) (b) 
Figure 10: The functions (a) X(x u l/R) and X (x t ) and (b) Y(x t ,l/R) and Y (x t ). 

200 GeV the functions X and Y are enhanced by 20% and 31%, respectively. For l/R = 250 GeV 
this enhancement decrease to 13%(21%) and are only 6%(9%) for l/R = 400 GeV. 

In table 2 we give the values of the functions S, C, X and Y for different l/R and m% = 
167 GeV. 



l/R [GeV] 


S 


C 


X 


Y 


200 


2.813 


1.099 


1.826 


1.281 


250 


2.664 


1.003 


1.731 


1.185 


300 


2.582 


0.946 


1.674 


1.128 


400 


2.500 


0.885 


1.613 


1.067 


SM 


2.398 


0.798 


1.526 


0.980 



Table 2: Values for the functions S 7 C, X and Y. 



4.8 Branching Ratios for Rare Decays 

The branching ratios for the rare decays in question can be directly obtained from [25] by simply 
replacing the SM functions Xq and Yq by X(x t , l/R) and Y{x t , l/R), respectively. We have 



Br(K 



7T ' vv) 



( ImAy 



V A 5 



X(x t ,l/R) 



+ 



H ' K MX) + ^X(x t ,l/R) 



X 



, (4.24) 



28 



3a 2 Br(K+ -> ir°e + v) „ ni 
«; + = r x+ % . . ^A 8 = 4.31 • 10~ n , (4.25) 



2vr 2 

where we have used [52] 



sin 



q = — , sin 2 9 W = 0.23, Er(K+ -► 7r°e + i/) = 4.87 • 10~ 2 . (4.26) 

Here Aj = V* s Vid with A c being real to a very high accuracy. = 0.901 summarizes isospin 
breaking corrections [53] in relating K + — > 7r + z/P to — > Tv e + v. -PoPO is given in (4.7). 
Next, 

Br(K L - tt ^) = KL • (^r*(zt, l/fl)) 2 , (4.27) 

KL = !^1^ K+ = 1.88. 10- 10 (4.28) 
r(K+) 

with K + given in (4.25) and rx L = 0.944 summarizing isospin breaking corrections in relating 
K L -► tt°i/P to if+ -► 7r°e+z/ [53]. 

Next, normalizing to Br(B — > X c eP) and summing over three neutrino flavours we find 

Br{B^2^1- 3a 2 |V^| 2 X 2 (x f , k(0) 

Br{B^X c eu) 4vr 2 sin 4 #J V cfe | 2 /(*) k(z) ' l ' j 

Here /(z) is the phase-space factor for B — > X c eP with z = m^/m 2 , and k(z) = 0.88 [54, 55] is 
the corresponding QCD correction. The factor k(0) = 0.83 represents the QCD correction to the 
matrix element of the b — > sz/P transition due to virtual and bremsstrahlung contributions. In 
the case of B — ► X^i/i/ one has to replace T4 S by which results in a decrease of the branching 
ratio by roughly an order of magnitude. In our numerical calculations we set f(z) = 0.54 and 
Br(B -► X c eu) = 0.104. 

Next, the branching ratio for B s — > is given by 



Br(g a ^Z + Z-)=r(^)^f - " ) F 2 m 2 m Bs Jl - 4^j-\V t * b V ts \ 2 Y 2 (x t , l/R) (4.30) 

7T \^ 4L7T Sill C'ly / y ^ 

where F^ s is the B s meson decay constant. The formula for Br{B& — > Z + Z~) is obtained by 
replacing s by d. The relevant input parameters are [32] 

F Bd = (203±§|) MeV, F Bs = (238 ± 31) MeV . (4.31) 

We set also r(B s ) = 1.46 ps and r(Bd) = 1.54 ps [52]. The short distance contribution to the 
dispersive part of — > [i + n~ is given by [46, 49] 

1 2 



-Br(K L -> /J~V )sd = 



^Po(F) + ^F(x t ,l/i?) 



(4.32) 
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= — — ^ A = L733 • 10 ' (433) 

where we have used Br(K + — ► = 0.634. The charm contribution including NLO correc- 

tions is given by [48] 

P (Y) = = 0.128 ±0.013. (4.34) 

Unfortunately due to long distance contributions to the dispersive part of Kl — ► n + /j>~, the 
extraction of Br(K^ — > ^ + /i _ )sD from the data is subject to considerable uncertainties. The 
present best estimate reads [56] 

5r(K L -► m + aOsd < 2.5 • 10- 9 . (4.35) 
4.9 GIM Mechanism and Convergence of the KK Sum 

The tree- level masses of the KK modes of all particles approach the value n/R for large KK 
mode numbers n, see eq. (2.30) and (A. 2). Due to the unitary CKM matrix the GIM mecha- 
nism [12] comes into play. It suppresses partly the higher KK mode contributions to the sums 
in eq. (3.4), (4.2), (4.5) and (4.10) and is essential for the determination of the dominant con- 
tributions. For the S n function the GIM suppression amounts to an additional factor of 1/n 4 
for the contributions WW, GG and WG. The contributions from diagrams with are not 
suppressed due to the second term in the coupling see (A. 3). This results in a hierarchy 
of the various contributions to S n with WW, GG and WG proportional to 1/n 6 , Wa and Ga 
proportional to 1/n 4 and the dominant contribution aa proportional to 1/n 2 for large values of 
n. 

For the penguin diagrams the effect of the GIM mechanism on the convergence of the partic- 
ular contributions is partly hidden by the subtle cancellation of divergencies among the different 
contributions and the two self-energy diagrams. For the combinations plotted in fig. 8 we ob- 
serve that the term corresponding to W ± : F\ + F2 + -F5 is logarithmically divergent without 
GIM and shows a 1/n 2 behaviour after GIM mechanism has been taken into account. The 
a ± , G^ : F3 + F4 + Fq + F-j term shows no GIM suppression and is proportional to 1/n 2 for 
large values of n. The mixed term W ± ,a ± ,G ± : Fg is constant for large n before the inclusion 
of GIM mechanism, but is GIM suppressed by a factor of 1/n 2 . 

The contributions to and B^ 1 show a 1/n 4 behaviour with and without GIM, except the 
Gww an d Hy/w terms. They are proportional to 1/n 2 without GIM and behave like 1/n 4 with 
GIM suppression at work. This different asymptotic behaviour with respect to the S n function 
is due to the appearance of leptons with negligible zero-mode masses instead of quarks in the 
box diagrams contribution to the functions G and H, see eq. (A. 3). 
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4.10 Numerical Analysis 

As discussed in section 3, \Vtd\, Q and f] in the SM and in the ACD model differ from each other. 
In our numerical analysis we will take this difference into account. Not taking this difference 
into account would misrepresent significantly the patterns of the enhancements of branching 
ratios in question. 




200 400 600 800 1000 200 400 600 800 1000 

R- 1 [GeV] R- 1 [GeV] 

Figure 11: Branching ratio for the decays K + — ► tt + vV and B s — > as predicted by the ACD model 

and the SM as functions of the inverse radius R of the extra dimension. 

In fig. 11 we show the branching ratios Br(K + — > n + uv) and Br{B s — ► as functions of 

the compactification scale 1/R. As IV^I^cd = l^tslsMj the enhancement of Br(B s — > fi + n~) is 
entirely governed by the ratio (Y/Yq) 2 . On the other hand the dependence of Br(K + — » -k + vv) 
on 1/R differs from the one of (X/Xq) 2 because of the additional charm contribution that is 
independent of 1/R and the fact that Vtd in the ACD model differs from its SM value. The 
remaining branching ratios are shown in table 3. The 1/R dependence of Br(B — > X s vv) is 
governed by the function X 2 (xt/1/ R), while the corresponding dependences of Br(B — > Xjyv), 
Br(K-L — > 7r°z/P), Br(Bd — > and Br(K^ — ► /i + /x~)sD include also the 1/R dependence 

of Vtd shown in fig. 4 (a). As expected, all the branching ratios are significantly enhanced for 
1/R < 300 GeV. For 1/R > 400 GeV, except for Br(B s -> the distinction between the 

predictions of the ACD model and the SM will be very difficult. 

In this numerical analysis we have used the results for the CKM parameters of fig. 4 and 
the central values of all the remaining input parameters as given above. The uncertainties in 
these parameters partly cover the differences between the ACD model and the SM model and 
it is essential to reduce these uncertainties considerably if one wants to see the effects of the 
KK modes in the branching ratios in question. Therefore a detailed analysis that includes all 
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l/R 


200 GeV 


250 GeV 


300 GeV 


400 GeV 


SM 


-r-» / 7- r 1 1 _\ -I /-> 1 1 

Br{K + — > 7r + ^z/) x lO 11 


8.70 


8.36 


8.13 


7.88 


7.49 


Br(K L — ► 7r u z/z/) x 10 11 


3.26 


3.17 


3.09 


2.98 


2.80 


Br(K\j -> )sd x 10 y 


1.10 


1.00 


0.95 


0.88 


0.79 


5r(B X^) x 10 5 


5.09 


4.56 


4.26 


3.95 


3.53 


Br(.B -► V d zvz>) x 10 6 


1.80 


1.70 


1.64 


1.58 


1.47 


5r(£ s -► /i+At~) x 10 9 


6.18 


5.28 


4.78 


4.27 


3.59 


Br(B d -► x 10 10 


1.56 


1.41 


1.32 


1.22 


1.07 



Table 3: Branching ratios for rare decays in the ACD model and the SM as discussed in the text. 



uncertainties would be in our opinion premature at present. 

4.11 An Upper Bound on Br(K + — > n + isis) in the ACD Model 

The enhancement of Br(K + — > ir + vi>) in the ACD model is interesting in view of the results 
from the AGS E787 collaboration at Brookhaven [57] that read 

Br(K + tt+vv) = (15.7^ 7 2 5 ) • lO" 11 (4.36) 

with the central value by a factor of 2 above the SM expectation. Even if the errors are sub- 
stantial and this result is compatible with the SM, the ACD model with a low compactification 
scale is closer to the data. As emphasized in [57, 58] the central value in (4.36) implies within 
the SM a value for \V t d\ that is substantially higher than the value obtained from the standard 
analysis of the UT of section 3. Here we would like to emphasize that within the ACD model 
Br{K + — > ir + vi>) is closer to the data in spite of the fact that [V^Iacd < IV^sm- The enhanced 
Z°-vertex represented by the function C is responsible for this behaviour. 

In [48] an upper bound on Br(K + — > tt + vu) has been derived within the SM. This bound 
depends only on \V c b\, Xq, £ and AM^/ AM S . With the precise value for the angle (5 now available 
this bound can be turned into a useful formula for Br(K + — > tt + ui>) [58] that expresses this 
branching ratio in terms of theoretically clean observables. In the ACD model this formula 
reads: 



Br(K + -► * + vv) = K + \V cb \ 4 X 2 (x t , l/R) 



aR 2 t sm 2 f3+-(R t cosP+ 



a V \V cb \ 2 X(x t , l/R) 

(4.37) 
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where a = 1/(1 — A 2 /2) 2 , R + = k + /X 8 and Rt is given in (3.29). This formula is theoretically- 
very clean and does not involve hadronic uncertainties except for £ in (3.29) and to a lesser 
extent in \V c b\. 

In order to find the upper bound on Br(K + — > 7r + i/p) in the ACD model we use 

\Vcb\ < 0.0422 , P (X) < 0.47 , sin/3 = 0.40 , m t < 172 GeV, (4.38) 

where we have set sin (3 to its central value (see (3.33)) as Br{K + — > ir + vv) depends very weakly 
on it. The bound on \V cb \ results from \V cb \ = 0.0406 ± 0.0008 [34]. We used here two standard 
deviations as the determination of \V c b\ involves some hadronic uncertainties. The result of this 
exercise is shown in table 4. We give there Br(K + — > Tr + uu) max as a function of £ and 1/R for 
two different values of AM S . The range for £ chosen by us is in accordance with (3.30) and the 
recent new analysis in [59] that gives £ = 1.22 ± 0.07. The upper bound in the SM given in the 
last column is lower than the values for l/R = 400 GeV by roughly 10%. We observe that for 
1/R = 200 GeV and £ = 1.30 the maximal value for Br(K + -> -k+vv) in the ACD model is 
rather close to the central value in (4.36). 

At first sight the 30% enhancement of Br(K + — > ir + vi>) for 1/R = 200 GeV with respect to 
the SM values seems to contradict the results in table 3, where a more modest enhancement of 
Br{K + — > tt + vu) is seen. However, one should realize that now AM<i/AM s and not AM^ alone 
enters the analysis and the enhancement of the function X is not accompanied by a suppression 
of \Vtd\i that with R t given by (3.29), equals the one in the SM. The consistency with the 
experimental value of AMj requires then a sufficiently small yj Bb a Fb a • In table 4 we indicate 

by a star the cases that require \J~B~b a Fb a < 190 MeV. Such low values are rather improbable 
in view of (3.30) and consequently values of Br(K + — ► ir + vi>) larger than 12 • 10~ n are rather 
unlikely even in the ACD model. 





1/R = 200 GeV 


1/R = 250 GeV 


l/i? = 300 GeV 


1/^ = 400 GeV 


SM 


1.30 


13.8* (12.3*) 


12.7* (11.3*) 


12.0* (10.7) 


11.3* (10.1) 


10.8 (9.3) 


1.25 


13.0* (11.6) 


12.0 (10.7) 


11.4 (10.2) 


10.7 (9.6) 


10.3 (8.8) 


1.20 


12.2* (10.9) 


11.3 (10.1) 


10.7 (9.6) 


10.1 (9.1) 


9.7 (8.4) 


1.15 


11.5 (10.3) 


10.6 (9.5) 


10.1 (9.0) 


9.5 (8.5) 


9.1 (7.9) 



Table 4: Upper bound on Br(K + — > ti + vv) in units of 10 11 for different values of £ and 1/R and 
AM S — 18/ps (21/ps). The stars indicate the results corresponding to J BB d Fs d < 190 MeV. 
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5 Summary and Outlook 



In this paper we have calculated for the first time the contributions of the Kaluza-Klein (KK) 
modes to AMk, £k, AM(j iS and rare decays K + — > tt + vv, — > 7r°z/P, i^L — ► ^ + ^ ^ B — > 
X s ^vD and S s ^ — > /x/Z in the Appelquist, Cheng and Dobrescu (ACD) model with one universal 
extra dimension. As a byproduct we have given a list of the required Feynman rules that have 
not been presented in the literature so far. 

The nice property of this extension of the SM is the presence of only a single new parameter, 
1/R. This economy in new parameters should be contrasted with supersymmetric theories and 
models with an extended Higgs sector. Taking 1/R = 200 GeV our findings are as follows: 

• The short distance one-loop function S(xt, relevant for AF = 2 transitions is larger 
than the corresponding SM function So(x t ) by roughly 17%. This implies on the basis 
of ek and AMd a 8% suppression of |T4d|ACD with respect to |Vtd|sM and a decrease of 
the angle 7acd by 10° with respect to 7sm- On the other hand, (AM s )acd is larger than 
(AM s )sm by 17%, see section 3 for details. AMk is essentially uneffected. 

• In order to see whether the modifications of the SM expectations are required by the data, 
the comparision of \Vtd\ and 7 extracted from sk and AMd in the SM and in the ACD 
model with the universal unitarity triangle constructed by means of \V u b/V c b\, AMd/ AM s 
and a(ipKs) will be important. To this end, the data on these quantities have to be 
improved and the uncertainties in the relevant non-perturbative parameters reduced. 

• The short distance one-loop function X(xt,l/R) relevant for the decays K + — ► tt + vu, 
Ki, — ► TT°uij and B — ► X 8 gyv is larger than the corresponding SM function Xo(xt) by 
roughly 20% due to KK contributions to the Z°-penguins. In the case of K + — > tt^vu, 
Kl — > 7r°vi> and B — ► X^vv this enhancement is partially compensated by the fact that 
l^tdUcD < |Vfci|sM and fj acd < VSM- We then find the enhancements of Br(K + — > ir + vv) : 
Br(K L -> 7r°i/p) and 5r(£ -> X d z/i>) over the SM expectations by 16%, 17% and 22%, 
respectively. As — ► X s vv is governed by the CKM element \Vt s \ that is common to the 
SM and the ACD model, the enhancement of Br(B — > X s vv) amounts to 44%. 

• The short distance one-loop function F(xt,l/i?) relevant for the decays Bd, s — ► //ju and 
ifx, — > /U/2, is larger than the corresponding SM function Yo(xt) by roughly 30%. As 
B s — > /i/Z is governed by the CKM element |Vj s | with (V^Iacd = I ^s|sm this implies a 72% 
enhancement of Br(B s — > //ju) relative to the SM expectation. In the case of 2?^ — ► /i/i and 
Kl — > /i/U, due to (V^Iacd < |^td|sMj the corresponding enhancements amount to 38% 
and 46%. 
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• As the ACD model belongs to the class of MFV models, general properties of these models 
identified in [14, 60] are automatically valid here. 

• For 1/R = 250 (300) GeV all these effects are decreased roughly by a factor of 1.5 (2.0). 
See fig. 11 and table 3. 

In short, the signatures of the deviations from the SM expectations are: 

• The decrease of \Vtd\, 7, fj and g. 

• The increase of AM S . 

• The increase of all branching ratios considered in this paper with a hierarchical structure 
of maximal enhancements: K + -> ir+vv (16%), K L -> tt'W (17%), B -> X d vv (22%), 
K L -► nfi (38%), B -► X s vv (44%), 5 d -► up (46%) and £ s -► nfi (72%) for 1/J2 = 
200 GeV. For l/R = 250 (300) GeV these enhancements are decreased roughly by a factor 
of 1.5 (2.0). 

In the coming years of particular interest will be the improved measurements of Br(K + — ► 
tt+vv) and AM S . Indeed, in the MFV models Br(K + — > 7r + i^) can be predicted as a function 
of V once the angle /3 and AMd/ AM S are known. The relevant formula is given in (4.37). For 
a given value of X, the branching ratio Br(K + — > 7r + z/p) decreases with increasing AM,. If 
the present central experimental values for Br(K + — > 7r + i/t>) will remain while the experimental 
error will decrease significantly, the SM expectations for Br(K + — »■ 7r + i/p) will be significantly 
below the experimental data while the corresponding estimates within the ACD model may 
agree with them due to X(x t , 1/R) > Xo(x t ), see table 4. 

Another definite prediction of the ACD model is an increase of AM S over the SM value. This 
should be contrasted with the prediction of the MSSM at large tan/3 where a suppression of 
AM S with respect to the SM value is predicted [61]. Thus, when the ratio (AM s ) cxp /(AM s )sm 
will be known with a sufficient accuracy, we will know whether the ACD model with a low 1/R 
or the MSSM with a large tan (3 is ruled out by the data. 

A distinction between the ACD model and the MSSM at low tan/3 will also be possible. 
In the latter case the super symmetric effects in AF = 2 transitions considered in section 3 
are generally larger than in AF = 1 transions so that the branching ratios for K + — > tt + uu, 
Kl — ► tt^uu, Kl — ► B — > X^i/v and -B^ — ► /i/Z are generally suppressed with respect to 

the SM while they can be enhanced for special ranges of supersymmetric parameters in the case 
of B —> X s vv and B s — > ///2 [62]. 

However, the main message from our analysis is the following one. Even for the lowest 
compactification scale, 1/R = 200 GeV, considered by us, the ACD model is consistent with 
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all the available data on FCNC processes analyzed here . No fine tunning of the parameters 
characteristic in the flavour sector for general super symmetric models is necessary. This is first 
of all connected with the GIM mechanism that assures the convergence of the sum over the KK 
modes in the case of Z° penguin diagrams, removing the sensitivity of the calculated branching 
ratios to the scale M s » 1/R at which the higher dimensional theory becomes non-perturbative 
and at which the towers of the KK particles must be cut off in an appropriate way. 

With the much improved data on the processes calculated in this paper and the theoretical 
uncertainties reduced, it should be possible in the future to distinguish the predictions of the 
ACD model from the SM ones, provided 1/R < 400 GeV. For higher compactification scales 
this distinction will be very difficult with a possible exception of B s — ► Whether these 

findings apply also to B — > X s j, B — > X s l + l~ and Kl — > ir°e + e~ is an interesting question. 
As the Z° -penguins are enhanced in the ACD model, the corresponding enhancement of the 
branching ratios for B — > X s l + l~ and Kl — > 7r°e + e~ can easily be calculated by means of 
the known formulae [25, 63, 64]. However, such an analysis would clearly be unsatisfactory as 
these decays receive also contributions from the 7-penguins and magnetic penguins. The KK 
contributions to 7-penguins are unknown, while the corresponding contributions to the magnetic 
penguins have been calculated in the context of an analysis of the decay B — > X s ^f under the 
assumption of the dominance of the scalars ajK [2] . We have seen that this assumption would 
be correct in the case of S(xt, 1/R) for 1/R > 300 GeV, but certainly would misrepresent the 
KK contributions to C(xt, 1/R) in view of strong cancellations between different contributions 
as discussed in section 4. Consequently a satisfactory analysis of B — > X s 7, B — > X s l + l~ and 
Kl — > 7r°e + e - requires a priori the inclusion of all KK contributions. We will address all these 
issues in the forthcoming paper [11]. 
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A Feynman Rules in the ACD Model 



In this section, we list all propagators and vertex rules needed for the calculation of the box and 
penguin diagrams considered in this paper. The Feynman rules are derived in the 5d i?g-gauge 
described in Section 2.3. 

It is convenient to define a 4 dimensional counterpart to every parameter in the 5 dimensional 
Lagrangian, marked there with a caret. The conversion factors are chosen in order to eliminate 



all explicit appearances of factors of \f2irR in the Feynman rules: 

fi = p,, v = V2kR v, g 2 = J— §2, 

1 ^ ^ (A.l) 

g, = vm Xu/V = Xu/V ' 

where p is the Higgs mass parameter. 

In order to simplify the notation, we omit the KK indices of the fields and of the mass 
parameters defined in (A. 3) and (A. 4). There is no ambiguity because in one- loop calculations 
at least one field is always a zero-mode. In our case, this is the Z boson in the Z vertices, and 
the down-type quark or the neutrino in all other vertices. Due to KK parity conservation, the 
other two fields have equal KK mode number, i.e. either zero or n > 1. 

Fermion zero-modes have substantially different Feynman rules than their KK excitations. 
The fields Q u , Qa, U, V, C u , C e and £ are always supposed to be (n > l)-modes, while the 
zero-modes are labeled u, d, v and e. The generation indices are i = u,c,t and j = d, s,b for the 
quarks and i,j = e,fi,T for the leptons. 

The masses of all bosonic particles can be expressed in terms of four independent mass 
parameters: 

t2 



M 2 



n 



(A.2) 

M l(n) = ^+V- 

We introduce two sets of mass parameters appearing in the fermion-scalar couplings. The 
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(A.3) 



mass parameters rrS^ are 
(i) n 

m l = ^ C i(n) +mSi( n) , 

(0 n 

™2 = -^Si(n) + r "i c i(n)> 

m 3 =-M w c i{n) + -—s i{n) , 

(i) n Trrii 

m\' = M w s i{n) + ^^; c i(n), 

where and the up-type fermion masses ra; on the r.h.s. are the zero-mode masses. The 
parameters Cj( n ) and Sj( n ) denote the cosine and sine respectively of the fermion mass mixing 
angle a^) defined in (2.29). The mass parameters Mx'^ are 

M i' j) = mjc i{n) , 



M 2 ,3) = m j s i(n), 

~~ c i(n) i 



(A.4) 



RM W 

Mr J> = ——J- 8i(n ), 

4 RM W l[n) ' 

where again My/ and the down-type fermion masses rrij on the r.h.s. are the zero-mode masses. 
All momenta and fields are assumed to be incoming. 



A.l Propagators 

The propagators are: 

for scalar fields S = H, a , a ± , A 5 , G°, G ± : 




k 2 — M 2 + ie 



with the masses 



s 


H(n) 


a (n) 






G \n) 


(n) 


M 


M H(n) 


M Z(n) 


M W(n) 









(A.5) 



for gauge bosons V = A,Z, W ± : 





k 




—i 








~ k 2 — M 2 + it 



<r - a - o 



k»k v 
k 2 - iM 2 + iej ' 
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with the masses 



V 



M 



A 



in) 



(n) 



w, 



(n) 



^A(n) ^Z(n) ^W(n) 



(A.6) 



for fermion fields F = u,d, Q u , Qd,U,V,u,e,C u , C e ,£: 

i($ + m) 
k 2 — m 2 + ie 

with the masses 




F 


u 


d 


Qu{n) 


Qd{n) 


U{n) 


V (n) 


V e C v (n) 




m 


m u 


m d 


m u(n) 




m u(n) 


m d(n) 




m e(n) m e(n) 



A. 2 Vertices 

The Feynman rules for the vertices are: 
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c w M W ( n ) 



ZW + G- : 


C = 


slM 2 ? 


9 n 
4- c — 
+ W R 2 ' 


ZW-G+ : 


C = 




2 n2 


ZW + a~ : 


C = 


„ , n 
~ W R 




ZW-a + : 


C = 


„ r n 





(A.7) 



(A.8) 

(A.9) 
(A.10) 
(A.ll) 
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ZG+G- 

Za + a~ 
ZG + a- 
ZG-a + 



(2 2 \ 2 2 

c w - s w) M W - 2c w jp, 



R 2 ' 



11 



C = M w -, 



C = -M w 



n 
R' 




i92C w (g^(k 2 - h) x + g^ x (h - k 3 ) u + g\ v (k 3 - k 2 




= |% (PlCl + PrCr) 



ZuiUi 



ZviUi : 



ZQiQi 



ZQMi : 



C L = 


3-44, 


Cr = 


-4s 2 


C l = 


3, 


Cr = 


o, 


Cl = 


-44 + 3c 2 () 


Cr = 


-44 + 3c 2 


C l = 


— 3Si(n) c j(n) j 


Cr = 


3Sj(n) c i(ra)i 



Zdjdj 



Zejej 



ZUiQi 



-3 + 2si 



Cr = 2sl 



C L 
Cr 

Cl 
Cr 

Cl 
Cr 



-3 + 64 



= 6* 



-44 + 3^( n )j 



-4s 2 + 3s 2 



(n)' 



= -3a 



i(n) Ci(n) 1 



= 3s 
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F 2 


F 1 

r 


W ± 
n 





l 92 „ r , 



W + Uidj : 


C L 


= Vij, 


w~ 


djUi : 


C L 


= v* 

v iji 


(A.21) 


W+ Qidj : 


C l 




w~ 


djQi : 


C l 


Cj(n) *ij i 


(A.22) 


W + Uidj : 


Cl 


= ~ s i(n)Vijj 


w~ 


djUi : 


Cl 


Tr* 


(A.23) 


W + V iej : 


C l 


= 


w- 


ejVi : 


Cl 


= Sij, 


(A.24) 


W + ViCj : 


C l 


= Sij, 


w~ 


CjVi : 


Cl 


= Sij, 


(A.25) 


W + Vi£j : 


C l 


= 0, 


w~ 


EjVi : 


Cl 


= 0. 


(A.26) 
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V2M, 



-(PlCl + PrCr). 



W(n) 



G + Uidj 
G + Qidj 



G + Uidj 



G + fiej 



G + ViCj 



G + Vi£j 



C L = -rriiVij, 
Cr = rrijVij, 

C L = -mfVij, 
C R = M^Vij, 

Cl = m^Vij, 
C R = -M?> j) Vij, 

C L = 0, 
Cr = rrijSij, 

C L = 0, 

(i) 

Cr = m\ J '5ij, 
C L = 0, 

Cr = -m^Sij, 



G djUi 



G- d,jQi 



G~ djUi 



G ejVi 



G Cj Vi 



G EjVi 



Cl = -mjV^ 
Cr = rmV^ 

C L = -Mf'%.*, 
C R = m?V* j , 

C L = M%' j) V* p 
C R = -mfv* p 

Cl = -rrijSij, 
Cr = 0, 

C l = -m^Sij, 
Cr = 0, 

C l = m^Sij, 
Cr = 0, 



(A.27) 
(A.28) 
(A.29) 
(A.30) 
(A.31) 
(A.32) 
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a + Qid 



a + Ujdj : 



a ViCj : 



a Vi£j 



C L 


W t r 




Cr 


= M^Vij, 


a 


C l 


(i) 

= myVij, 




Cr 


= -M^Vij, 


a~ 


C l 


= o, 




Cr 


= 777,3 8ij, 


a~ 


C l 


= o, 




Cr 




a~ 



"3 ^ 



"3"i 







Cr = 


rr$V*, 


C l = 




Cr = 




C L = 




Cr = 


o, 


Cl = 




Cr = 


0. 



(A.33) 

(A.34) 
(A.35) 

(A.36) 



B Different Contributions to AF = 2 Box Diagrams 

The contributions from different sets of diagrams to the functions F(x t ( n ), x u ( n )) in (3.8) are 
given as follows 



F WW(n) = U { x t(n), x u(n)), 

W(n) 



(B.l) 



n ) 

W(n) 

x ^(aJt(n),aJ u (n)), 



(*) , (*) l r («) i w 

m l c t(n) ' m 2 s t(n) m l c u(n) ~r m 2 s u(n) 



(B.2) 



n 

*Wa(n) = -2 Tjg 

1V1 W(n) 

x U(x t ( n ), X M ( n )), 



(*) , (*) l r («) , («) 

™3 C <(n) + m A s t(n) m\ C M ( n ) + m\ s u{n) 



(B.3) 



1 



Ga(n) 2 ^ 6 



(*) (*) , (*) (*) 



l W(n) 



(u) (u) (u) (u) 

m\ 777,3 + m 2 m 4 



F GG(n) 



1 MS 



4M« x 

W{n) 



{mff + {mff (mn a + «) a C/(^ (n) , x u(n) ), 



(B.4) 
(B.5) 



W(n) 



{mf? + {mf? KY + KT tf(s,(„),s«<„)). 



( M )\2 1 /V„( u )\2 



(B.6) 
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The functions U and U are denned as follows 



U(x t ,x u ) 



x\ logxt 



+ 



x 2 log x u 



+ 



(x t - x u )(l - x t ) 2 (x u - x t )(l - x u ) 2 (1 - x u )(l - x t y 



2x t \ogx t , l + x t 
U{x t ,x t ) = — — + 



il-x t f {l-x t f 



U(x t ,Xi 



x t log x t 



+ 



X u log X u 



+ 



(x t - x u )(l - x t ) 2 (x u - x t )(l - x u ) 2 (l-x u )(l-x t ) 



\ (1 + at) log gt , 
U(x t ,x t ) = — - - V 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



(l-x t ) 3 (l-^) 2 ' 

C Different Contributions to Z°-Penguin Diagrams 

The contributions of the diagrams in fig. 7 to the functions F{xf) in (4.12) are given as follows 

F l(Xf(n)) 



1 / 4 , 4 4 2 

C /(n) + S /(n) - 3*u- 



A + 111 + K ( X f(n)) ~\- 2Xf(n)hq (x /(n) ) 



^2(x /(n )) 
^/(n)) 



1 2 2 
4 C /W S /(n) 



A + In ■ 



16m 



W(ra) 



- - + h q (x /(n) ) + 2z /(n) /i, (x /(n) ) 

m /(n) ^ 

rff + (-F) 2 ) (<$<„) " |-«) + ((4 /} ) 2 + (mi 



A + 111 ^ + ^9 fc/w) - 2Xf(n)hq (s/(„)) 



^4(^/(n)J - -^T72 l m l m 2 + m 3 m A ) C f(n) s f{n) 

SIV1 W(n) V ' 



A + ln n+ h i ( x f(n)) + 2x f (n)h q (x f(n) ) 

m f(n) Z 



F 5(Xf(n)) 

F&(x f{n) ) 



3 -c 2 
4 Cw 



A + ln 



/i- 



1 



ML, , 6 

Vv (n) 



16m W(n) 



(l-2 4)M^ + 2c^) (rfV + rf) 2 ) 
+ ((1 " 2 4) ^ + 24M 2 ,) ((m^) 2 + (mf ) 2 )" 

A + ln j£ + \~ X f(n)K (x /(n) ) 



(C.l) 

(C.2) 



2 _ 41 2 
s /(n) 3 S «; 



(C.3) 



(C.4) 
(C.5) 



(C.6) 
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F 7(Xf(n)) 



Here 



8M4 R V'" 1 

W(n) 



A + ln + \~ X f(n) hw (*/(„)) 

VT(n) 



*s(a;/(n)) = 



/(»») 



Vv (n) 



2 , ,r2 2 n2 \ ( (f) (f) \ 

s w M w -c w -j^J [m\' c /(n) + m\ s f(n) ) 



+M]V ~R~ C/(n) +m 4 /} S /(n)) ^ ( X /(n)) • 



(C.7) 



(C.8) 



A = - + ln4vr - 7£, D = 4 - 2e 
e 

and the functions /i g and h w are given by: 



h q (x) 
h w (x) 



1 



lnx 



1 — x (1 — x) 2 ' 

1 xlnx 
+ 



(C.9) 

(CIO) 
(C.ll) 



l-x ' (1-x) 2 ' 

Finally, the contributions from counter terms corresponding to the self-energy diagrams of fig. 
12 that should be added to the functions are given by 



1 / l + x f{n) 2x 2 m \nx m 



X f(n) (1 " X f{n) f 



In 



M W(n) 



AS2 {Xfin)) = 8 ^ \ 2 I " (1 _ x/(n)) 2 



In- 



M ¥ 2 



W(n) 



(C.12) 
(C.13) 




G ± ,a d 



j Qi, Ui 



(1) 



(2) 



Figure 12: Self-energy diagrams neccessary for calculating the clcctroweak counter term as discussed in 
[50]. 
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D Closed form for the AC(x t ,l/R) function 



We express the logarithms of (4.14) as integrals I\ and I2 

x„ln^±^ = h(x t )-h(l)+x t -l, (D.l) 
In^i = h{x t )-I 2 {l) (D.2) 

1 T X n 

with 

h(a) = a 2 C dy 2 V , (D.3) 
Jo ay + x n 

h(a) = a [ dy . (D.4) 

Jo ay + x n 

Next, we interchange integration and summation. The integrands can now be summed using 
the relation 

Eb _b (^/c7r coth(^/c7r) - 1) 

n=l 

This allows us to derive a closed form for the sum 



E° (x t - 7) x t _ mw^Rxt 

, C * x » x »> ~ 16(1-*,) 16(1 -x t f 



n=l 

x I dy 
Jo 



3(l + x t ) 



coth(mw^R\/y) — xl coth(m t 7ri?y / y) I 



+ (x t — 7)y/y (coth(mw^R^/y) - xl coth(mtirRy/y) 



(D.6) 

Expanding in [mw^Ry/y^t^Ry/y) and integrating (D.6) leads to the expression in (4.17). 

E Different Contributions to AF = 1 Box Diagrams 

The one loop amplitude of the KK excitations in the diagrams in fig. 9 is a sum of contributions 
coming from the various bosonic fields in the loop 

G n {Xf{n)i x e{n)) = G\VW{n) + G WG ( n ) + GVa(n) + G Ga („) + G GG r( n ) + G aa ( n ), (E.l) 
H n (Xf(n)> x v(n)) = HwW{n) + %G(n) + %a(n) + ^Ga(n) + ^GG(n) + ^aa(n)- ( E - 2 ) 
The explicit results for B^ l/ (xt,x n ) and Bn^(x t ,x n ) are given by 

-B n (xt,x n ) = G n (x^( n ), x e ( n )) — G n (x u ( n ), x e ( n )), 

n)i x u(n)) ■^■ni x u(n) i x u(n 

)), (E.3) 
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where the results for the vv-box diagrams take the form 

Gw\V(n) 
G\VG(n) 



—2 U(x f(n) ,X e(n) ), 

1V1 W{n) 



W m f(n) m e(n) 



G 



Wa(n) 



2 M* f , 

W(n) 

l M W m f(n) m e(n) 

W(n) 



m 



i )c f(n) + ™-2 >s f(n) rnf U{x f(n) ,x e{n) ) 



,(/) 



m i f)c f(n) + m ( fsf(n) ™ff U{x f(n) ,x e{n) ), 



G 



Ga(n) 



G 



GG(n) 



G, 



aa(n) 



1 Myy 

W(n) 

1 Myy 
16M.6 

W(n) 

1 Mw 

W(n) 



(f) (f) (f) (f) 



mf^mj f7(x /(n) ,x e(n) ), 



(m^) 2 + (m^) 2 ] (mS e) ) 2 C/(x /(n) , x e(n) ), 



(m^) 2 + (m^f] (^ e) ) 2 f/(x /(n) ,x e(n) ), 



and for the /j/i-box diagrams we get 



H 



WW(n) 



1 M 2 , 



4M2 . 

W (n) 



#1 



VKG(n) 



ff«(n) 



1 M> /( 

" ) n 

2 Mf, /( " 

Vv (n) 

1 Jl^m^m,^) 
2 



m[ f) c f{n) +m { 2 f) s f{n) mf ] U{x f(n) ,x u{n) ), 



Vv (n) 



m 3 c /(«) + m 4 ' s /(n) m !T ^/(n) , «i/(n) )> 



,(/) 



(") 



Ga(n) 



1 M 2 , 

W(n) 
1 Ml 



(t) (f) (f) (f) 

m\ 777.3 + m 2 m X 
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i ; m 3 C/(x/(„),x„( n )), 



GG{n) 



W 



Vv (n) 



(777S /) ) 2 + (m( /) ) 2 l rf } ) 2 C/(x /(n) ,x, (n) ), 



if, 



1 M 2 



aa(n) 



W 



16M.6 

Vv (n) 



(m^) 2 + (mf) 2 ] (m^) 2 t/> /(n) , *„(„)). 



(E.4 
(E.5 
(E.6 
(E.7 
(E.8 
(E.9 

(E.10 
(E.ll 
(E.12 
(E.13 
(E.14 
(E.15 



The functions U and U are defined in (B.7)-(B.10). We have taken into account the overall 
minus sign in (4.8). 

Summing all the contributions we find 



p 17x t + 18x n x t - 9x n x n (26x t + 9x n + 18x n x t ) x n 
B n (x t ,x n ) = — 1 In 



16(x t -l) 16xt 
(x n + x t ) (9x n + 17xt) j X n + Xt 
16(x t - lfx t 1 + x n ' 



1 + X, 



(E.16) 
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3x n + 5xt Qx n xt x n (3x n + 6x n x t - 2x t ) 



16(x t -l) ' ' \Qx t 

- hx t ) (x n + x t ) , x n + x t 
H n In 

16(x t - lfx t 1 + X n 



1 + X n 



(E.17) 
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